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Background 


9 NASA has identified cryogenic storage & transfer as an 
area with greatest potential for cost saving 

9 Zero-Boil-Off (ZBO)or Reduced Boil-Off (RBO) dynamic 
storage tank pressure control that involves some mode 
of mixing of the bulk liquid with or without active or 
passive cooling is needed to realize cost savings. 

9 Optimization of ZBO or RBF tank design for 
microgravity applications will most probably be 
accomplished with only ground-testing due to 
budgetary constraint. 

9 State-of-the-Art validated storage tank CFD models will 
play an crucial role in extrapolation of the lg tested 
storage tank design to microgravity and partial gravity 
applications. 



9 Correct implementation of Interfacial & bulk turbulence and evaporative condensing mass 
transfer is crucial for the fidelity and validity of the CFD models. 
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Fundamental Multiphase Science Issues 





Natural Connect (turbulence) 
Forced Mixing (turbulence) 
Evaporation Condensation 
Microgravity Superheats 


Non-Condensable Gases 
Transport Barrier 
Marangoni Convection 
Interfacial Kinetics 


Free Surface Dynamics 
Contact Angle Dynamics 
Sloshing/Droplet Transport 
Phase Control/Positioning 




Two-Phase Sharp Interface Storage Tank CFD Model 
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Continuity: 


^ + vOcv) = 0 
dt ^ ’ 

Momentum: 

— (pv)+V(pvv) = -Vp + V[p € (Vv + Vv r )]+ 


Interfacial Energy Balance: 

LJ V = -fej V7 ; • n + k-VT-n V 



Schrage Interfacial Mass Transfer: 



Energy: 
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Two-Phase VOF Storage Tank CFD Model 

Continuity: ^ + v(pP) = 0 

dt 

& 

Momentum: — (pv)+v(pvv) = -Vp + V \Meff ( V V + Vl ’ T )]+ PS + F, 



|-( pE)+ v(y(pF + p))=v(k eff vr)+ s „ 


Energy: _d_ 

Volume of Fluid (VOF) model: 

Energy and Temperature are defined as mass average scalars: 
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Properties: 
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Continuity of Volume Fraction of the q-th phase 
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Interfacial mass transfer per unit volume: 

s a = ril, • A. 

m -sec 

A . = Vcr , is an interfacial area density in 1/m, 1X1, is a mass flux vector in kg/(m 2 -sec). 

where a is a volume fraction of the primary phase 

Schrage's Relation : 
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MHTB Facility (1996-2005) 



9 Self-pressurization tests were run with LH2 

9 Cylindrical midsection with: 

height = 3.05 m 
diameter = 3.05 m 

9 Ullage pressure, 2 Fluid temperature rakes, 
Large number of Wall Temperature 
measurements 

9 Boil-off test was performed prior to tank 
lockup and self-pressurization 

9 Most tests include 20, 50, 90% fill levels 

9 Heat Flux = 2.05 W/m 2 

9 Tank Internal volume 37.5 m 3 

9 2:1 elliptical top and bottom domes 

9 Tank is enclosed in a vacuum shroud. 


Vacuum Chamber 


\ 



Multi-purpose Hydrogen Test Bed (MHTB) 
at NASA Marshall Space Flight Center 
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Effect of Accomodation Coefficient 


MHTB Pressurization - MHTB 90% 
Sharp Interface Model 




Heat Flux = 2.05 W/m 2 


MHTB Sharp Interface 90% fill ratio, Turbulent, AC=1 
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P, Pa 



Effect of Accommodation Coefficient - MHTB 90% 

Sharp Interface Model 



MHTB 90% Self-Press with Sharp Interface - Schrage model 

Heat Flux = 2.05 W/m 2 
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Effect of Accommodation Coefficient- (MHTB 90%) 

Sharp Interface Model 
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Heat Flux = 2.05 W/m 2 
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Pressure, Pa 



Effect of Accommodation Coefficient - MHTB (50%) 

VOF Model 




time after beginning of self-pressurization, seconds 
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Vapor Heat Flux = 0.90 W/m 2 
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Mass Transfer Rate at the interface, kg/s 



1G MHTB Pressure Control Results: LH2, Large Tank 

Droplet Spray Bar - VOF 
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Effect of Accommodation Coefficient - Summary 



9 Computation of mass transfer using Shrage in the sharp Interface 
model is insensitive to the magnitude of accommodation coefficient, 
a close to 1 seems to work fine. 

9 Computations of mass transfer using Schrage in the VOF method also 
show a certain degree of insensitivity to the magnitude of 
accommodation coefficient but small values around a = .01 seem to 
be necessary for numerical practicality and stability 

9 This is only true for a stable flat interface between bulk phases. Mass 
transfer computations based on Schrage equation for droplets and 
during slosh dynamics or boiling situations are quite sensitive to 
magnitude of accommodation coefficient. 

9 Schrage might not represent the right mass transfer kinetics under 
these conditions. 
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Effect of Turbulence 


Effect of Turbulence: MHTB (50%) 




Liquid Heat Flux = 2.05 W/m 2 
Vapor Heat Flux = 0.90 W/m 2 
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Pressure, Pa 


MHTB Self Pressurization: Effect of Turbulence 


Pressure Time History 
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MHTB Self-Pressurization: Effect of Turbulence 
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Tank Self-Pressurization Experiments at K-site Facility 

( 1990 - 91 ) 



1 . Test fluid is liquid hydrogen 

2. Flightweight insulated 2219-T62 aluminum ellipsoidal tank 

• Internal volume: 4.95 m 3 = 175 ft 3 

• Tests conducted in vacuum chamber. 

• Test article is enclosed by a cryoshroud whose 
temperatures are maintained with electrical heaters. 

• Tank is insulated with 2 blankets of MU. 

3. Steady boil-off test and measurement performed at 95% 
liquid fill fraction and 117 kPa (or 1.17 bar) tank pressure. 

4. Tank fill level was reduced to desired fill level (29%, 49%, 
83%) 

5. Several hours of additional venting at 103 kPa were 
performed to achieve stationary state. 

6. Self-pressurization tests were initiated from a stationary 
stratified state. 

7. Two Cryoshroud Temps -> Two heat loads (2 & 3.5 W/m 2 ) 

8. Grashof Number (Gr) based on 3.5 W/m 2 average heat flux 
into tank -> vapor: Gr = 2.21e+13; liquid: Gr = 1.33e+14 
(which corresponds to turbulent natural convection for a 
steady-state natural convection flow) 
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P, Pa 


Effect of Turbulence: K-Site (50%) 


K-Site VOF, Conjugate HT, Laminar, 3.5 W/m A 2 


Tank Pressure & Interfacial Mass 
Transfer 




T[K] 

43.10 

37.68 

32.27 

26.85 

21.43 



21.43 



32.27 


26.85 


27.40 


21.43 1 


22.04 


K-Site VOF, Conjugate HT, Turbulent, 3.5 W/m A 2 
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60,000 seconds of self-pressurization 


K-Site Self-Pressurization: Effect of Turbulence 


Results: 3.5 W/m 2 heat flux 
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K-Site Self-Pressurization: Effect of Turbulence 




Results: 3.5 W/m 2 heat flux 
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K-Site Self-Pressurization: Effect of Turbulence 


Results: 2.0 W/m 2 heat flux 


Tank Pressure 


I Experiment 
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K-Site Self-Pressurization: Effect of Turbulence 


Velocity Magnitude 
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Effect of Turbulence- Summary 




9 Laminar models agree closely with the pressure evolution and vapor 
phase temperature stratification but under-predict liquid 
temperatures. 

9 Turbulent SST k-w and k-e models under-predict the pressurization 
rate and extent of stratification in the vapor but represent liquid 
temperature distributions fairly well. 

9 These conclusions seem to equally apply to large cryogenic tank 
simulations as well as small scale simulant fluid pressurization cases. 

9 Appropriate turbulent models that represent both interfacial and bulk 
vapor phase turbulence with greater fidelity are needed. 

9 Application of LES models to the tank pressurization problem can 
serve as a starting point. 
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K-Site Self-Pressurization: Effect of Turbulence 

Sharp Interface Model 
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